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APERTURE  AVERAGING  OF  SCINTILLATION 
FOR  SPACE-TO-GROUND  OPTICAL  COMMUNICATION  APPLICATIONS 


Communication  from  space-to-ground  stations  which  utilize  coherent  short 

wavelength  radiation  (~lym)  as  the  carrier  must  contend  with  atmospherically 

induced  scintillation.  The  amount  of  scintillation  with  which  the 

0 

communication  system  must  deal  is. a  function  of  the  diameter  of  the  collecting 
entrance  pupil  in  the  ground  station  optical  system.  As  is  well  known. 


increasing  the  size  of  the  aperture  diameter  reduces  the  amount  of  scintilla¬ 


tion  measured  by  the  detector 


1,2,3 


The  aperture  averaging  factor  is  defined  as  the  ratio  of  the  variance  of 
irradiance  obtained  from  a  finite  size  collecting  aperture  to  the  correspond¬ 
ing  quantity  obtained  from  a  "point  aperture."  In  this  letter  the  aperture 

averaging  factor  is  calculated  for  both  spherical  and  plane  waves.  The 

2 

Kolmogorov  spectrum  with  an  arbitrary  index  structure  constant  Cn  profile,  as 

would  be  encountered  in  a  space- to-ground  propagation  geometry,  is  used 

here.  A  useful  engineering  formula  for  the  aperture  averaging  factor  that 

exhibits  the  explicit  parametric  dependence  on  wavelength,  collector  diameter 

2 

and  path  weighted  integrals  of  is  presented  below.  In  Ref.  3  it  is  shown 
that  the  aperture  averaging  factor,  A,  can  be  written  as 


.  1  C,(Dx) 

A  ■  s  /  "error  Vx)  *** 

o  I 


where  C^(p)  is  the  covariance  of  irradiance  and 


M,  (x)  ■  (2/ir){cos  *x  -  x  (1-x2)  ] 


is  the  optical  transfer  function  of  an  unaberrated  circular  entrance  aper¬ 
ture.  Here  we  consider  weak  scintillation  conditions  only  where  Cj(0)  «  1. 

Following  the  development  given  in  Ref.  3,  it  can  be  shown  for  spherical 
waves  that  the  aperture  averaging  factor  can  be  written  as 


cx(o) 


where 


Z  oo 

B  -  4s2  /  dz'(z-z*)2£-)2  /  f  f_*  K5dK 


z  '  '  s  R  n 

o 


oo 

Cx(0)  -  4ir2  /  dz’<z-z’>2£-)2  /  fs*nK5dK 


f  »  sine2 [K2z' (z-z* )/2kz] 
s 


2Jx(KDz ’/2z) 


and,  for  the  Kolmogorov  spectrum, 


♦j^K,z ' )  *  0.033  C2(z')K“11/3, 
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In  Eqs.  (3)-(6),  z'  ■  0  corresponds  to  the  location  of  the  source,  z  is  the 
distance  between  the  source  and  the  receiver,  k  is  the  optical  wavenumber, 
sinc(u)  *  sin(u)/u,  D  is  the  diameter  of  the  circular  collector,  and  is  the 

Bessel  function  of  the  first  kind  of  order  one.  The  quantities  fg  and  fR  can 
be  regarded  as  the  spherical  wave  turbulent  diffraction  and  receiver  filter 
functions,  respectively.  To  modify  Eqs.  (3)  -  (6)  for  plane  wave  propagation 
the  factor  z'/z  is  replaced  by  unity. 

We  are  primarily  concerned  with  downward  propagation  from  an  exoatmos- 
pheric  laser  source  to  a  receiver  located  within  the  atmosphere  or  on  the 
ground.  For  the  corresponding  case  of  upward  propagation,  the  transverse 
irradiance  correlation  scale  length  is  typically  much  greater  than  the  dia¬ 
meter  of  collection  apertures  of  practical  concern,  and  hence  no  appreciable 
aperture  averaging  effects  will  result  (i.e.,  A  *  1). 

For  spherical-wave  propagation,  fR  possesses  spatial  frequencies 
K  <  KR  -  2z/Dz '  a  2/D.  The  relative  importance  of  turbulent  diffraction 
effects  and  aperture  averaging  can  be  seen  by  comparing  Kr  and  Kg,  where  Kg  « 

[2kz/z'(z-z')]1/2  »  [2k/z  ]l,l  For 

o 

D2  »  ^  ,  (8) 

Kr  «  Kg  and  all  the  turbulence  lies  in  the  near  field  of  the  aperture.  In 

this  case  the  geometrical  optics  formulation  (i.e.,  f  “1)  can  be  used  in 

s 

Eq.  (3)  to  compute  the  collected  scintillation.  On  the  other  hand,  for  D2  << 
zQ/k  it  follows  that  Kr  »  Kg,  and  one  may  use  the  approximation  fR  <*  1  in 
Eq.  (3)  with  the  result  that  the  aperture  averaging  factor  approaches  unity. 
The  quantity  zQ  is  of  the  order  10  km  and  thus,  for  example,  apertures  of  dia¬ 
meters  larger  than  a  few  inches  will  exhibit  averaging  effects  for  visible 
light  propagation. 


In  the  downward  propagation  direction,  the  atmospheric  turbulence  is 

limited  to  the  latter  portion  of  the  propagation  path  where  z'  *  z.  The 
2 

function  C^Cz')  is  nonzero  only  in  the  range  z  -  zq  <  z'  _<  z,  where  for  an 

exoatmospheric  source  zq  «  z.  A  change  of  variables  from  z'  to 

2  2 

n  *  z  -  z*  is  indicated  with  C^z')  replaced  by  Cn(n)  *  where  n  *  0  corre¬ 
sponds  to  the  location  of  the  receiver.  If  inequality  (8)  is  satisfied, 

then  f  *  1  in  Eq.  (3)  and  the  aperture  averaging  factor  is  given  by 
s 

,  7/6 

A  a  A  (Ah  lu  )  (9) 

o  o 


where 


A  =  4 
o 


r  -2/3  .  2,  .  , 

J  x  Jj  (x)  dx 

l _ 

r  -H/16  ,  2  . 

J  x  sin  x  dx 


r5(i/3)  _ 

14 >3  1/2  ~  0.90, 

nWJ[31/ Ml 


(10) 


h 
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path 


/  dnC2(n)n2( 1-nz  1) 


-1/3 


,16/7 


/  dnC^(n)[n(l-nz  *)] 


path 
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A  is  the  optical  wavelength,  and  r  is  the  gamma  function.  The  quantity  hQ 
can  be  regarded  as  an  atmospheric  turbulence  aperture  averaging  scale 
height.  For  plane  waves  or  for  spherical  waves  where  z  »  zQ 

(i.e.,  n/z  «  1)  we  have  that 


h 

o 


path 


/  dn  C2(n)n2 


/  dn  C2(n)n5^6 


6/7 


(12) 


An  important  attribute  of  Eq.  (9)  is  that  it  reveals  the  explicit  parametric 

dependence  of  the  aperture  averaging  factor  on  optical  wavelength,  collector 

2 

diameter,  and  path  weighted  integrals  of  C  .  For  propagation  at  zenith 

n 

angle  0  it  is  easy  to  show  that  hQ  in  Eq.  (9)  is  replaced  by 

2 

hQsec  0,  and  Cn(n)  in  Eqs.  (11)  and  (12)  is  now  given  by  the  index  structure 
constant  profile  for  the  situation  of  interest. 

An  engineering  formula  can  be  constructed  from  Eq.  (9)  which  can  be 

useful  even  when  inequality  (8)  is  violated  (i.e.,  for  small  apertures). 

Indeed,  for  «  zQ/k,  A  «  1,  and  thus  without  requiring  detailed  knowledge 

of  Cj(p),  the  aperture  averaging  factor  can  be  approximated  by  the  engineering 
formula: 

A.  - — 1 - m  <13, 

1  +  A  [D  /Ah  sec  0] 
o  o 

where  A  ^  =  1.1.  This  expression  gives  the  aperture  averaging  factor  for 

space-to-ground  propagation  conditions  in  closed  form  and  is  thus  suitable  for 

system  studies  of  a  broad  class  of  atmospheric  turbulence  models. 

For  illustrative  purposes  we  present,  in  Table  1,  numerical  values  of  hQ 

o 

and  A  for  zenith  propagation,  A  =  1  pm,  D  ■  1  m,  and  various  profiles 
quoted  in  the  literature.  In  addition,  we  give  the  corresponding  weak  scin¬ 
tillation  values  for  the  "point  aperture"  variance  of  irradiance: 
cx(o)  »  2.24  k7/6/cJ(n)n5/6dn. 

Examination  of  Table  1  reveals  for  the  turbulence  models  considered  here 
that  both  h^  and  A  based  on  the  NAVY/DARPA  model  are  about  a  factor  of  two  to 


four  less  than  the  corresponding  quantities  based  on  Barleti's  and  Hufnagel's 
models.  For  example,  the  NAVY/DARPA  daytime  model  and  the  Hufnagel  model  (for 
V  *  27  m/sec)  both  indicate  that  the  variance  of  irradiance  at  a  "point"  is 
about  0.1.  On  the  other  hand,  the  ratio  of  the  corresponding  aperture  averag¬ 
ing  factors  is  about  3.6.  Hence,  for  X  *  1  pm  and  aim  entrance  aperture 
diameter,  the  variance  of  irradiance  based  on  the  Hufnagel  model  is  3.6  times 

larger  than  that  obtained  from  the  NAVY/DARPA  daytime  model.  As  this  example 

2 

illustrates,  care  should  be  exercised  in  the  application  of  available  Cp 
models  to  optical  communication  system  performance  studies. 


Turbulence 

Model 

Scale  Height 
hQ  (km) 

Aperture  Averaging  Factor: 
Zenith  Propagation, 

X  =■  1pm,  D  *  lm 

Point 

Irradiance 

Variance 

navy/darpa(4) 

(DAYTIME) 

3.4 

1.2  x  10-3 

0.095 

NAVY/DARPA 

(NIGHTIME) 

5.3 

2.0  x  10“3 

0.053 

Barletl  et  al.^3^ 

9.8 

4.1  x  10"3 

0.18 

hufnagel(6) 

(V  ■  18  m/sec) 

9.1 

3.7  x  10“3 

0.063 

HUFNAGEL 

(V  ■  27  m/sec) 

10.3 

4.3  x  10"3 

0.11 

HUFNAGEL 

10.8 

4.6  x  10"3 

0.19 

HUFNAGEL 

(V  ■  36  m/sec) 


10.8 


0.19 
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LABORATORY  OPERATIONS 


The  Laboratory  Operations  of  The  Aerospace  Corporation  la  conducting  exper¬ 
imental  and  thaoratlcal  investigations  nacaaaary  for  tha  evaluation  and  applica¬ 
tion  of  aclantlflc  advances  to  oew  military  apaca  systems.  Versatility  and 
flexibility  have  bean  developed  to  a  high  degree  by  the  laboratory  peraonnel  in 
dealing  with  the  aany  problems  encountered  In  the  nation's  rapidly  developing 
•pace  systems-  Expertise  In  the  latest  scientific  developments  la  vital  to  the 
accomplishment  of  tasks  related  to  these  problems.  The  laboratories  that  con¬ 
tribute  to  this  research  are: 

Aerophyslcs  Laboratory:  Launch  vehicle  and  reentry  aerodynamics  and  heat 
transfer,  propulsion  chemistry  and  fluid  mechanics,  structural  mechanics,  flight 
dynamics;  hlgh-temperature  thermomechanics,  gas  kinetics  and  radiation;  research 
In  environmental  chemistry  and  contamination;  cw  and  pulsed  chemical  laser 
development  Including  chemical  kinetics,  spectroscopy,  optical  resonators  and 
beam  pointing,  atmospheric  propagation,  laser  effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions,  atmo- 
apherlc  optics,  light  scattering,  state-specific  chemical  reactions  and  radia¬ 
tion  transport  in  rocket  plumes,  applied  laser  spectroscopy,  laser  chemistry, 
battery  electrochemistry,  space  vacuum  and  radiation  effects  on  materials,  lu¬ 
brication  and  surface  phenomena,  thermionic  emission,  photosensitive  materials 
and  detectors,  atomic  frequency  standards,  and  bloenvlronmental  research  and 
monitoring. 

Electronics  Research  Laboratory:  Microelectronics,  CaAs  low-noise  and 
power  devices,  semiconductor  lasers,  electromagnetic  and  optical  proportion 
phenomena,  quantum  electronics,  laser  communications,  Ildar,  and  electro-optics; 
communicstlon  sciences,  applied  electronics,  semiconductor  crystal  and  device 
physics,  radiometric  imaging;  mllllmeter-vave  and  microwave  technology. 

Information  Sciences  Research  Office;  Program  verification,  program  trans¬ 
lation^  performance-sens 1 t lve  system  design,  distributed  architectures  for 
spaceborne  computers,  fault-tolerant  computer  systems,  artificial  intelligence, 
and  microelectronics  applications. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metal  matrix 
composites,  polymers,  and  new  forms  of  carbon;  component  failure  analysis  and 
reliability;  fracture  mechanics  and  stress  corrosion;  evaluation  of  materials  in 
space  environment;  materials  performance  in  space  transportation  systems;  anal¬ 
ysis  of  systems  vulnerability  and  survivability  in  enemy-induced  environments. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radiation 
from  the  atmosphere,  density and  composition  of  the  upper  atmosphere,  aurorae 
and  alrglov;  magnetospherlc  physics,  cosmic  rays,  generation  and  propagation  of 
pl***a  waves  in  the  magnetosphere;  solar  physics.  Infrared  astronomy;  the 
effects  of  nuclear  explosions,  magnetic  storms,  and  solar  activity  on  the 
earth's  atmosphere,  ionosphere,  and  magnetosphere;  the  effects  of  optical, 
electromagnetic,  and  particulate  radiations  in  space  on  space  systems. 


